AD-A03J  124 


UNCLASSIFIED 


TRACOR  INC  AUSTIN  TEX  -/a  1T/1 

£5  2iST2RJx2Hc£«A  JTL^2LJrLCCTC0  FR0"  THE  0CEAN  **™~etJ<u> 

MAR  65  H P BUCKER*  ® I BOuRX ANOFF  N0002^*67"C«1S1A 

TRACOR -6A- 34 3~<J  C ‘i*8 


IDA  033 1 24 


10^1  81  — (fjU  ■ 


Contract  N00024-67-C-1318 
Project  Serial  SF101-03-16 
Task  8524 

TRACOR  Project  002  046  01 
Document  Number  68-343-U 


*.  most  Project 

/ft 


TECHNICAL  MEMORANDUM 


THE  DISTORTION  OF  A SONAR  PULSE 
REFLECTED  FROM  THE  OCEAN  SURFACE 


H.  P.  Bucker 


G.  I.  Bourianoff 


Submitted  to: 


Commander , Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Attn:  Code  00V1C 


March  1,  1968 


DISTRIBUTION  STATEMENT  A 

Approved  ioi  public  Twhtnut 
Pitributivo  Unlimited 


6500  Tracor  Lane,  Austin,  Texas  78721,  AC  512/926-2800 


t 


ntmniiiii/irJm 


NOV  26  1976 


^vU  U 


W TRACOR 


A 


810303  0071 


Thai:  oh 


6S00  TRACOR  LANE.  AUSTIN.  TEXAS  78721 


ContractnNflQ024-67-C-131fl 
Project  Serial  SF101-03-16 
Task  8524 

TRACOR  Project  002  046  01 
Document  Number  68-343-U 


TECHNICAL  ME! 


THE  DISTORTION  OF  A SONAR  PULSE 
REFLECTED  FROM  THE  oHeAN  SURFACE 


Bucker 


Bourianoff 


Submitted  to 


Commander,  Naval  Ship  Systems  Conmand 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Attn:  Code  00V1C 


March  1,  19( 


Approved 


Submitted: 

U-  P. 

H.  P.  Bucker 
Project  Director 


W.  C.  Moyer 
Program  Manager 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  78721 

ABSTRACT 

Equations  are  derived  for  the  direct  evaluation  of  a 
sonar  pulse  reflected  from  a time-varying  ocean  surface.  Fewer 
assumptions  and  approximations  are  made  in  this  direct  approach 
than  in  conventional  developments  of  scattering  theory.  However, 
the  cost  of  computer  time  may  become  excessive  for  some  cases. 

The  calculated  surface  reflected  signals  are  correlated 
with  the  incident  sonar  pulse  to  find  the  expected  correlation 
loss  of  a sonar  processor.  For  a short  (5  msec)  low  frequency 
FM  slide  the  calculated  correlation  losses  (1.6  dB  for  the  sur- 
face generated  by  a 10  knot  wind  and  3.9  dB  for  a 20  knot  wind) 
were  almost  entirely  the  result  of  energy  loss  in  the  reflected 
signal. 
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1.  INTRODUCTION 

In  most  operations,  an  active  sonar  system  must  process 
a signal  which  has  reflected,  in  part  or  whole,  from  the  ocean 
surface.  The  effect  of  signal  distortion  because  of  the  surface 
reflection  should  be  understood  if  the  processing  system  is  to 
do  a good  job  of  identifying  possible  target  echoes. 

Many  investigators^"”'*  have  made  experimental  measure- 
ments of  backscattering  coefficients.  Forward  scattering  coeffi- 
cients are  more  difficult  to  measure  because  the  source  and 
receiver  must  be  separated.  Measurements  have  been  reported  by 
Liebermann^  and  Marsh^. 

Many  theoretical  studies  have  been  reported  where  equa- 
tions are  derived  relating  the  scattering  coefficient  to  sea 

7-12 

state  parameters  for  a stationary  surface  . All  of  these 
theories  assume  the  incident  wave  to  be  single  frequency,  con- 
tinuous, and  plane.  Also,  all  of  these  theories  contain  approxi- 
mations that  become  invalid  as  the  surfaces  become  rough  or  the 
grazing  angles  become  small.  As  would  be  expected,  the  theoreti- 
cal predictions  are  accurate  for  large  grazing  angles  and  reason- 
ably smooth  surfaces  but  have  limited  usefulness  under  other 
conditions. 

Model  scattering  experiments  have  been  reported  by 
La  Casce  and  Tamarkin*"'*,  by  Horton,  Mitchell,  and  Barnard^, 
and  by  Medwin*^.  These  have  been  most  helpful  in  testing  the 
different  theoretical  approaches  to  the  scattering  problem.  Their 
usefulness  in  helping  to  understand  the  active  sonar  problems  has 
been  limited,  however,  because  of  the  difficulty  in  scaling  the 
model  results  to  the  full  size  ocean  environment. 


The  reflection  of  a sinusoidal  wave  from  a time  varying 
surface  has  been  studied  by  Liebermann^,  Marsh ^ and  Parkins  ®. 
Liebermann  demonstrated  the  presence  of  Doppler  shifted  components 
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in  the  reflected  signal.  Marsh  made  modifications  to  his  earlier 
work  so  that  the  Doppler  shifts  could  be  calculated.  Parkins 
developed  the  theory  in  such  a manner  that  explicit  forms  of  the 
spectral  density  of  the  reflected  waves  were  obtained  in  the 
limiting  cases  of  slightly  rough  and  very  rough  surfaces. 

The  above  references  are  only  a small  fraction  of  the 
papers  concerning  scattering  theory.  They  were  listed  because 
they  represent  significant  advances  in  the  field  and  because  they 
indicate  the  current  points  of  attack  on  the  problem.  However, 
in  none  of  the  above  references,  or  in  any  other  literature 
reviewed,  is  any  method  for  accurately  calculating  the  distortion 
of  a sonar  pulse  wh^n  it  reflects  from  this  time  varying  ocean 
surface.  To  answer  this  question,  work  was  begun  at  TRACOR  to 
develop  a direct  calculation  of  the  reflected  sonar  pulse  using 
Kirchoff's  integral.  Details  of  the  calculation  are  given  in 
Section  4. 

Although  the  assumptions  usually  made  in  scattering 
theory  were  not  required  in  the  equations  developed  here,  the 
present  method  is  limited  by  its  running  time  on  a computer. 

The  running  time  is  proportional  to  the  number  of  reflecting 
areas  on  the  surface  and  to  the  length  of  the  pulse.  In  order 
to  study  long  pulses  (a  500  nsec  3.5  kHz  pulse  contains  over 
1700  cycles),  the  surface  insonification  has  been  limited  in  the 
present  work  to  a strip  one  meter  wide  with  the  strip  length 
determined  by  the  3 dB  down  points  of  the  beam  pattern.  It  seems 
plausible  that  this  narrow  beam  will  have  about  the  same  distor- 
tion as  a more  typical  sonar  beam,  but  this  supposition  remains 
to  be  proved. 

Section  3 describes  the  generation  of  the  ocean  sur- 
face. The  surface  is  formed  by  a set  of  100  moving  cosinusoidal 
waves  with  random  phases.  Choosing  a new  set  of  phases  will 
result  in  a completely  different  surface,  but  one  that  has  the 
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same  energy  spectrum.  Thus  the  two  surfaces  are  different  samples 
from  the  same  sea. 

In  Sections  5 and  6 some  typical  reflected  signals  are 
displayed  and  the  distortion  in  the  reflected  pulse  is  determined 
by  means  of  replica  correlation.  Correlation  losses  that  would 
occur  in  typical  sonar  processors  are  calculated. 

Section  7 describes  an  attempt  to  improve  the  operation 
of  a correlator  by  modifying  the  replica  to  fit  the  distortion  in 
the  reflected  signal.  Although  the  improvement  gained  in  the 
correlator  output  was  nil,  this  result  has  significance. 

Section  8 describes  a simple  method  of  setting  a general 
type  of  pulse  (such  as  a pseudo  random  noise  pulse)  in  a form 
suitable  for  efficient  evaluation  of  Kirchoff's  integral. 
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2 . SUMMARY 

Equations  have  been  derived  and  a computer  program 
written  for  the  calculation  of  a sonar  pulse  that  reflects  from 
the  time-varying  ocean  surface.  No  restrictions  on  the  size  of 
the  waves  or  type  of  sonar  pulse  are  inherent  in  the  program, 
although  the  cost  of  computer  time  may  become  excessive  for  some 
cases . 

To  calculate  the  sonar  distortion  for  a typical  ping, 
the  pulse  form,  beam  shape,  and  wind  speed  are  specified.  The 
program  then  generates  a time  varying  sample  of  the  ocean  sur- 
face and  the  reflected  pulse  is  calculated  by  a direct  evalua- 
tion of  Kirchoff's  integral.  The  program  can  be  run  again  with 
the  same  input  parameters  but  with  a different  sample  of  the 
surface.  Thus,  ping  to  ping  variation  of  the  reflected  pulses 
can  be  studied  in  a quantitative  manner. 

The  effect  of  the  surface  reflection  is  separated  into 
two  parts.  First,  the  fluctuation  of  energy  in  the  reflected 
pulse  for  successive  pings,  i.e.,  for  different  surface  realiza- 
tions, is  tabulated.  Second,  the  amount  of  correlation  loss 
because  of  signal  distortion  is  calculated  by  comparing  the 
correlator  output  for  the  reflected  signal  with  the  correlator 
output  for  the  pulse  reflected  from  a flat  surface.  For  a short, 
5 msec,  FM  slide  the  correlation  loss  for  a 10  knot  wind  speed 
was  1.6  dB,  all  of  which  was  due  to  energy  loss.  The  correlation 
loss  for  a 20  knot  wind  speed  was  3.9  dB,  of  which  3.8  dB  was 
due  to  energy  loss. 

Equations  are  developed  for  matching  the  correlator 
replica  to  the  distorted  pulse,  however,  the  correlator  improve- 
ment was  negligible.  Finally,  a generalization  is  discussed 
that  will  permit  calculation  of  reflected  pseudo  random  noise 
pulses  and  reflected  explosive  pulses. 
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3.  SIMULATION  OF  THE  SEA  SURFACE 

The  moving  sea  surface  is  simulated  by  superimposing 

a large  number  of  cosinusoidal  waves  of  random  phase.  This 

model  is  called  a long-crested  Gaussian  linear  Eulerian  model  of 
19 

random  waves  . The  elevation  of  the  sea  surface,  z,  as  a func- 
tion of  distance,  in  the  direction  the  waves  are  traveling,  and 
time  is 

z(x,t)  » V a2n+l  cos (k2n+lx_u,2n+lt  + e 2n+l^  ' 
n-1 

2 

For  deep  water  k is  equal  to  uu  /g,  w is  the  angular  frequency, 

g is  the  gravitational  constant,  and  the  {e}  are  a set  of  random 

phases  chosen  from  a distribution  uniform  over  0 to  2tt,  The 

coefficients  of  the  traveling  cosine  waves,  a,  are  related  to 

19 

the  spectral  density  of  the  sea  surface  S(uu)  by 

a2n+l  " [2S(u)2n+l)  ' (u,2n+2  " ®2n^  * 

Using  semi -empirical  arguments  as  to  the  form  of  S (ou)  and  a large 

20 

number  of  experimental  data  points  Pierson  and  Moskowitz  derived 
the  following  formula  for  S(u>)  : 

S(iu)  - (ag2/u>5)  exp  (-B  (g/(vuu))4)  . 

o 

Here  u>  is  the  angular  frequency,  g is  9.8  m/sec  , v is  the  wind 
velocity  in  m/sec,  and  a and  e are  dimensionless  constants  with 
a - 0.0081  and  0 - 0.74.  Figure  3-1  is  a plot  of  S(oj)  vs  fre- 
quency for  some  typical  values  of  wind  speed  in  knots. 

Samples  of  the  sea  surface  at  different  times  are 
generated  by  superimposing  100  traveling  cosine  waves  with  random 
phases  shown  in  Figs.  3-2,  -3,  -4,  -5.  Figures  3-2,  3-3,  and 
3-5  were  calculated  using  the  same  set  of  random  phases  but  cor- 
respond to  wind  speeds  of  5,  10,  and  20  knots  respectively.  The 
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four  curves  in  each  of  the  figures  represent  the  surface  at  four 
consecutive  times  separated  by  \ sec  intervals.  Figure  3-4  cor- 
responds to  a 10  knot  wind  speed  but  for  a different  set  of  ran- 
dom phases.  Thus  Figs.  3-3  and  3-4  correspond  to  different 
samples,  or  representations,  of  the  same  sea. 
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4.  MATHEMATICAL  FORM  OF  THE  REFLECTED  PULSE 

In  this  section  an  equation  is  derived  for  calculation 
of  the  reflected  pulse  as  a function  of  the  incident  pulse  form, 
the  geometry  of  the  problem,  and  the  motion  of  the  surface.  Re- 
ferring to  Fig.  4-1,  the  source  at  x , z emits  a pulse  described 

s s 

in  terms  of  a velocity  potential  f(t)D(0).  The  pulse  font  is 
determined  by  f(t)  while  D(0)  takes  care  of  the  beam  pattern  and 
level  of  the  source.  The  equation  of  the  moving  surface  is 
z = z(x,t) . 

Consider  the  problem  if  the  surface  were  removed.  The 
potential  at  any  point  p would  be 

• =f(ct  - p)  • D(0)/p  , (4.1) 

where  p is  the  distance  from  the  source  to  p and  c is  the  velocity 

of  sound.  With  the  surface  in  place  the  potential  at  p,  any 

point  in  the  water  or  at  the  surface  can  be  written  as 

f(t  - p /c)  • D(0 ) /p  + U.  Because  U is  a potential  due  to  the 

presence  of  the  surface  it  may  be  called  the  scattered  potential 

function.  Since  the  acoustical  pressure  must  be  zero  at  the 

surface, 

f(ct  - p)  • D(0)/p  + U = 0 (at  the  surface)  . (4.2) 

Q 

Following  Eckhart  , the  normal  derivative  of  U,  dU/dn, 
is  equal  to  plus  or  minus  the  normal  derivative  of  the  incident 
field  depending  upon  whether  the  paths  from  the  source  and 
receiver  to  the  surface  point  are  clear  or  shadowed.  Thus 


" q °(0)  &(f(ct  - p))  , (4.3) 

where  q is  +1  for  a clear  surface  area  and  -1  for  a shadowed  area. 
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FIG.  4-1B 


FIG.  4-1  - GEOMETRY  OF  SCATTERING  SYSTEM 
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The  complete  and  rigorous  solution  for  the  reflected 

21 

sound  field  can  be  written  in  terms  of  the  Kirchoff  integral 


IT  1 r f r ||  1 S / 1\  1 5 rr  ^ 3-1  1 r ~ S-\  1 

Up  " Sn^P  ” cr  d"ir^dt  ^ “ r ^n  ^ 


au. 


dS. 


(4.4) 


Here  U is  the  velocity  potential  of  the  scattered  pulse  at  a 
point  P,  and  U is  the  potential  at  the  surface.  The  integral 
is  sunmed  over  the  insonified  surface  and  the  square  brackets  [] 
mean  that  the  potential  functions  should  be  retarded.  Referring 
to  Eqs.  (4.2),  (4.3)  and  Fig.  4-1  it  is  not  difficult  to  show 
that  at  the  surface, 


Ug  « - f(ct  - p)/p 

[U8]  - - f (ct  - s)/p  , s - r + p 

dU_/dn  - q d { f (ct  - p)/p}/dn 
8 

5/Sn  •-  - cos  a S /dp 

a Ug/d t - -(i/p)  a f(ct  - p)/at  - (-c/p)f' 

3 r/d  n - - cos  p 

d(l/r)/ar  - cos  p/r2  (4.5) 

Substituting  (4.5)  into  (4.4),  and  with  some  more  algebraic 
reduction,  Kirchoff' s integral  can  be  written  as 


0 (t)  - - 2£1  r (t  g (£°SjL  + 

■ Q — * “ 


pr 


P r 


+ [f'](£2£-6.  + 1 ££?  g))  dS  . 


(4.6) 
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The  significance  of  the  two  values  of  q can  be  made 
clear  by  using  Eq.  (4.6)  to  calculate  the  pulse  reflected  from 
a perfectly  flat  pressure  release  surface.  It  can  be  shown 
that  if  q is  equal  to  one, the  reflected  pulse  is  a perfect  replica 
of  the  incident  pulse  but  with  change  of  sign,  which  is  the  same 
answer  as  from  the  method  of  images.  If  q is  -1  the  reflected 
pulse  is  identically  zero  which  is  the  appropriate  return  from 
a shadowed  section  of  the  rough  surface.  Finding  the  shadowed 
areas  is  a straight-forward  problem;  however,  in  the  cases  con- 
sidered so  far  all  surface  elements  have  been  assumed  to  be 

clear.  An  estimate  of  the  fraction  of  scattered  areas  can  be 

22 

made  using  equations  derived  by  Wagner 

Any  attempt  at  direct  evaluation  of  Eq.  (4.6)  leads  to 
considerable  difficulty.  For  example,  for  a given  t the  location 
and  orientation  of  each  surface  element  must  be  found  at  the 
earlier  time  at  which  the  surface  element  reflected  the  sound 
energy  which  arrives  at  the  receiver  P at  t.  A much  easier 
approach  is  to  focus  attention  on  a single  surface  element,  for 
example,  the  surface  between  x^  and  x ^ in  Fig.  4-2,  and  calculate 
the  complete  reflected  signal  from  this  region  as  a function  of 
time.  The  region  of  interest  is  then  shifted  to  an  adjacent 
area,  the  reflected  signal  calculated  and,  using  the  principle 
of  superposition,  the  second  signal  is  weighted  according  to  Eq. 
(4.6)  and  added  to  that  already  calculated.  When  all  regions 
have  been  covered  the  complete  reflected  signal  has  been  calculited. 
A simple  example  will  illustrate  the  details  of  the  calculation. 

Consider  a pulse  made  up  of  two  one-cycle  sine  waves, 
with  periods  T^  and  T^,  as  shown  in  Fig.  4-2*.  The  time  when  the 

*tn  this  section  we  will  assume  that  f(t)  can  be  represented  by  a 
set  of  sinusoidal  cycles.  The  period  of  each  cycle  may  be  chosen 
as  any  function  of  time.  This  provides  an  excellent  representation 
of  CW  pulses  and  FM  slides.  In  a later  section  it  will  be  shown 
how  any  reasonable  form  of  f(t) , including  pseudo-random  noise 
pulses  and  explosive  pulses,  can  be  built  up  of  segments  of  sinu- 
soidal waves  plus  a DC  component.  The  computer  time  using  the  more 
general  pulse  will  not  be  appreciably  greater  than  for  the  pulses 
considered  here. 
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REGION  I 


REGION  2 


rf(ct-r)  * sin[k,(ct - r ) J , k,  = 2ir/(cT,) 
M'(ct  -r)  » k,cos(ct-  r) 

rt  ( c t - r)  * *in[k2(  c(  t - T, ) - r )]  , k2  * 2r/(c  T2) 
V'(ct  - r ) s k2  cos[c(t  - Tj)  - r] 


FIG.  4-2  * SAMPLE  CASE  FOR  TWO  CYCLE  PULSE 
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beginning  of  the  first  cycle  reflects  from  region  1 is  found  as 
follows.  If  the  surface  were  flat,  this  time  would  be  t°  * p°/c, 
where  o°  is  the  distance  between  the  source  and  a flat  surface 
in  region  1.  Let  z°  be  the  actual  surface  elevation  at  this  time 
and  z°  be  the  z component  of  surface  velocity.  Then  the  correc- 
tion to  the  time  of  reflection  is  |z  |z°/(o°c  - |z  | z°) , where  z 

s s s 

is  the  z coordinate  of  the  source.  The  time  of  arrival  of  the 
beginning  of  the  first  cycle  at  the  receiver  is  p°/c  + r°/c  + 6t^, 
where  r°  is  the  distance  between  the  receiver  and  the  flat  surface 
in  region  1 and 

6tl  “ I 2SI  Z°/<P °c  " l2s|z°)  + |zr|z°/(r°c  - | zr|  z°)  . 

The  same  procedure  is  repeated  to  find  the  beginning  of  the 
second  cycle  (time  - p°/c  + r°/c  + T^  + 6t2)  and  the  end  of  the 
second  cycle.  Now  with  the  values  of  cos  a,  cos  0,  and  dS,  for 
the  times  when  the  sound  waves  reflected  from  the  surface  element, 
the  complete  wave  reflected  from  region  1 is  determined.  The 
wave  is  sampled  and  the  sampled  values  are  stored  in  a string  of 
computer  cells  according  to  the  sampling  time.  The  whole  procedure 
is  repeated  for  the  other  surface  regions  and,  using  the  principle 
of  superposition,  the  sampled  values  of  the  signal  from  each  new 
region  are  put  in  on  top  of  previously  calculated  values  in  the 
computer. 
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5.  CALCULATED  SURFACE  REFLECTED  PULSES 

The  final  objective  of  this  project  is  to  calculate 
the  correlation  loss  and  energy  fluctuations  that  an  active 
sonar  will  encounter  because  of  surface  reflections.  However, 
the  computer  runs  to  date,  which  will  be  reported  here,  have 
been  short,  5 msec,  FM  slides  with  frequencies  around  1 kHz. 

The  short  pulse  length  obviously  conserves  computer  time  and 
the  low  frequency  allows  a visual  inspection  of  the  distortion 
in  the  reflected  pulse.  A sampling  rate  of  10  kHz  is  built  into 
the  program  so  that  sonar  systems  operating  at  frequencies  up 
to  5 kHz  can  be  studied. 

The  most  interesting  result  from  the  calculations  to 
date  is  that  the  replica  correlator  does  an  excellent  job  of 
identifying  the  distorted  signal.  Although  ping  to  ping  fluctua- 
tions of  energy  levels  in  the  reflected  signals  are  large  and 
the  reflected  signals  show  obvious  distortion,  the  correlator 
output  has  no  loss  because  of  distortion  in  the  reflected  signal 
for  a 10  knot  wind  speed  and  only  0.1  dB  loss  for  a 20  knot 
wind . 

Figure  5-1  shows  the  incident  pulse  used  in  this  study 
and  the  pulse  reflected  from  a flat  ocean  surface.  The  incident 
pulse  is  a Doppler  invariant  (i.e.,  the  period  of  each  cycle  is 
a linear  function  of  time)  FM  slide  with  5 cycles  sweeping  from 
950  to  1050  Hz.  The  source  is  located  50  feet  below  the  surface 
and  sends  out  a 10  degree  vertical  beam  at  a grazing  angle  of  20 
degrees.  The  receiver  is  also  50  feet  below  the  surface  and 
positioned  in  the  line  of  the  specular  reflection. 

Figure  5-2  shows  the  incident  pulse  and  five  reflected 
pulses  calculated  by  using  five  different  representations  of  a 
sea  generated  by  a 10  knot  wind.  For  the  computer  print-out  the 
largest  sample  of  each  reflected  wave  was  normalized  to  one. 
Compared  to  a total  energy  of  1.0  for  the  pulse  reflected  from 
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FIG.  5-1  TRANSMITTED  AND  REFLECTED  PULSE 
FROM  A FLAT  OCEAN  SURFACE 
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a flat  surface  the  relative  energies  of  the  five  reflected 
signals  are  0.77,  0.99,  0.78,  0.79,  and  0.55.  In  the  next  sec- 
tion the  values  are  converted  to  dB  levels  with  the  mean  and 
standard  deviation  listed  in  Table  6-1. 

Figure  5-3  shows  the  incident  and  reflected  pulses  for 
a 20  knot  wind  speed.  The  relative  energy  ratios  for  the  reflected 
waves  are  0.29,  1.24,  0.13,  0.58,  and  0.09. 
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6.  DISTORTION  IN  THE  REFLECTED  PULSE 

After  the  reflected  signal  is  calculated  the  question 
arises,  by  what  criterion  should  be  distortion  be  measured? 

Here,  distortion  is  defined  in  terms  of  the  loss  of  correlation 
that  an  optimum  linear  processor  would  suffer  because  of  the 
distorted  form  of  the  reflected  signal. 

It  is  well  known  that  the  optimum  linear  processor  for 

identifying  a known  signal  in  white  noise  (i.e.,  produces  maximum 

S/N  ratio  at  output)  can  be  realized  as  a matched  filter  or  as 

23 

a linear  replica  correlator  . In  the  present  work  the  Fast 
Fourier  Transform  method  is  used  to  calculate  the  correlation 
between  the  surface  reflected  signal  and  the  incident  signal. 

Let  {yk,  k = 0,  1,  2,  ...,  K-l]  represent  K equally 
spaced  samples  of  the  signal  reflected  from  the  surface  and 
fx^,  k = 0,  1,  2,  ...  K'-l)  represent  samples  from  the  incident 
signal.  Since  the  reflected  signal  is  extended  in  time  by  the 
reflection  process,  K is  greater  than  K' . The  correlation  of 
the  two  signals,  Cg,  corresponding  to  a time  lag  sT,  T = time 
interval  between  samples,  is  given  by 

K-JL-s 

cs  - r *k  ^ • <6a> 

k-0 

For  long  pulses,  the  time  required  for  evaluating  Eq.  (6.1)  can 
be  substantially  reduced  by  using  the  Fast  Fourier  Transform. 

The  details  of  the  Fast  Fourier  method  will  be  discussed  next 
to  prove  that  the  correct  answer  is  obtained  and  to  illustrate 
the  required  arrangement  of  { x^]  and  { y^} . 

Let  be  a sequence  of  2N  numbers,  the  first  K'  of 
which  are  the  same  as  x^  with  the  remaining  terms  set  equal  to 
zero.  The  value  of  N is  specified  by  requiring  that  N * 2P, 
where  p is  the  smallest  integer  such  that  2P  * K.  Let  be 
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a sequence  of  2N  numbers,  the  first  N of  which  are  zero,  the 
next  K of  which  satisfy  the  relation  YN+k  = yk,  with  the  last 
remaining  terms  set  equal  to  zero.  Let  A^_  be  components  of  the 
discrete  Fourier  transform  (DFT)  of  (Xk)  and  Br  be  components 
of  the  DFT  of  fYj  . That  is, 

2N-1 

a = Y wrk 

Ar  L Xk  W 
k=0 


B 

r 


2N-1 

L 

t-0 


where , 


W s exp(-  tt i/N) 


(6.2) 


Next,  form  the  sequence  {D^]  where  Dr  is  the  product  of  complex 
conjugate  of  times  B^.  The  inverse  discrete  Fourier  trans- 
form of  f Dr)  can  then  be  written  as 

2N-1 


1 

2CT 


B W 
r 


-rs 


1 r r r 

m Li  l 

r k l 


X._  Y.  W 


-rk  ,.,r<t 
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■rs 


1 


k t 


Xk  Y. 


wr(t-k-s) 

/ t 

r 


2N-1-S  2N-1 

Z Xk  Vs  + Z Xk  Yk+s-2N 

k-0  k- 2N-s 


(6.3) 


2N-1 

The  last  line  follows  from  the  fact  that  W**1  equals  0 when 

r-0 

n is  an  integer,  except  when  n - 0,  ±2N,  *4N,  ...  in  which  case 
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the  sum  is  equal  to  2N.  The  last  sum  is  zero,  because  of  the 

zeroes  in  {X^l  and  fY^l . The  other  sum  in  the  last  line  is  the 

desired  correlation  function  C . The  above  arguments  are  similar 

s24 

to  those  used  by  Cochran  et  al  . in  the  evaluation  of  convolu- 
tion integrals. 

The  auto-correlation  function  for  the  incident  signal 
is  shown  at  the  top  of  Fig.  6-1  and  the  five  cross-correlation 
functions  for  the  reflected  pulse  at  10  knot  wind  speed  are  shown 
below.  All  of  the  correlograms  are  normalized  to  a peak  value 
of  one.  The  correlator  output  is  squared  and  averaged  over  a 
5 msec  period  in  a detector.  The  detector  outputs  are  shown  in 
Fig.  6-2.  The  correlograms  corresponding  to  a 20  knot  wind  speed 
are  shown  in  Fig.  6-3.  The  detector  outputs  are  similar  in  form 
to  those  in  Fig.  6-2  and  are  not  shown. 

The  effect  of  the  surface  reflection  on  processor 
performance  is  summarized  in  Table  6-1.  First  there  is  a per- 
formance loss  because  the  energy  level  is  down.  This  is  shown 
in  column  1^  where  the  reference  energy  level  is  that  of  the 
signal  reflected  from  the  flat  surface.  The  next  column  is  the 
standard  deviation  of  the  energy  levels  divided  by  the  five- 
sample  mean  energy  level. 

The  total  correlation  loss  is  shown  in  column  J5.  This 
was  calculated  by  comparing  the  average  correlator  output  energy 
for  the  signal  reflected  from  the  time-varying  surface  to  that 
for  the  signal  reflected  from  the  flat  surface.  The  fluctuations 
in  this  correlation  loss  are  given  in  column  6. 

The  total  correlation  loss  can  be  separated  into  two 
parts.  Part  is  due  to  the  energy  loss  and  part  results  from  the 
distortion  of  the  signal.  The  correlation  loss  due  to  distortion 
is  the  difference  between  the  total  loss  and  the  energy  loss. 

This  is  shown  in  column  3_.  As  shown,  the  total  loss  for  a 10 
knot  wind  is  1.6  dB  which  is  all  due  to  energy  loss.  At  20 
knot  wind  speed  the  total  loss  is  3.9  dB,  3.8  dB  because  of 
energy  loss  and  0.1  dB  because  of  distortion. 
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FIC.  6-1  - CROSSCORRELATION  FUNCTIONS  OF  THE  SIGNALS 
SHOWN  IN  FIG.  5 2 CORRELATED  WITH  THE  TRANS- 
MITTED SIGNAL  (10  KNOT  WIND  SPEED) 
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FIG. 6-2  - CROSSCORRELATION  FUNCTIONS  OF  FIG.  6*1 
AFTER  SQUARE  LAW  RECTIFICATION  AND  PER- 
FECT AVERAGING  (10  KNOT  WIND  SPEED) 

27  ^QjjjJ^owc  **-140- 14 

2 /?<✓«•  •oJwAmorr-tj 


STANDARD  DEVIATION 
OF  TOTAL  CORRELATION 
LOSS  DIVIDED  BY  MEAN 


MEAN  TOTAL  COR- 
RELATION LOSS 


STANDARD  DEVIATION 
OF  DISTORTION  COR- 
RELATION LOSS  DI- 
VIDED BY  MEAN 


MEAN  CORRELATION 
LOSS  DUE  TO  SIGNAL 
DISTORTION 


STANDARD  DEVIATION 
OF  ENERGY  DIVIDED 
BY  MEAN 


MEAN  ENERGY  LEVEL 


00 

.296 

GO 

GO 

O 

Q 

vO 

<7\ 

• 

• 

■ 

GO 

eg 

LTV 

GO 

O 

CT\ 

O 

O 

GO 

CD 

O 

O 

O 

— 

• 

• 

bJ 

<_> 

UJ 

Ui 

< 

o 

o 

ti- 

< 

< 

er 

Li. 

ti- 

=> 

er 

er 

CO 

=5 

3 

CO 

CO 

t- 

o 

t- 

»- 

z 

o 

r 

z 

-J 

r 

u. 

o 

o 

•OUftlANOFF 


6S00  TRACOR  LANE  AUSTIN.  TEXAS  78721 


7.  DISTORTED  REPLICA  CORRELATION 

Because  small  improvements  in  the  figure  of  merit  of 
a sonar  system  will  generally  result  in  a large  increase  of  the 
maximum  detection  range,  great  effort  is  spent  to  see  that  each 
part  of  the  system  operates  with  maximum  efficiency.  However, 
there  are  fundamental  limits  on  each  component.  For  example, 
the  onset  of  cavitation  restricts  the  maximum  transducer  radiation 
power  while  the  physical  size  of  the  transducer  limits  the  direc- 
tivity factor  of  the  beam.  For  the  signal  processing  sections 
of  the  sonar  system  the  theoretical  limits  of  performance  are  not 
as  well  set. 

In  this  report,  the  signal  processing  consists  of  cor- 
relation plus  detection.  The  correlator  multiplies  the  received 
signal  y(t)  by  a time  shifted  replica  signal  r(t-r).  Here  t is 
the  time  shift  or  time  lag  and  r(t)  is  generally  set  equal  to 
the  signal  sent  out,  denoted  x(t) . It  is  well  known  that  this 
type  of  correlator,  called  the  linear  replica  correlator,  is  the 
optimum  (i.e.,  produces  the  largest  signal-to-noise  ratio)  linear 
system  if  the  returning  signal,  y(t),  has  the  form  y(t)  z.  const. 

• x(t)  + white  noise.  In  this  section  a method  is  developed  for 
adjusting  the  replica  to  fit  the  distortion  of  the  returning 
signal,  thus  improving  the  correlator.  Other  possibilities  for 
modifying  the  replica  would  include  increasing  the  pulse  length 
by  adding  new  cycles  or  adding  low  frequency  components  which 
would  alter  the  envelope  of  the  pulse.  In  Fig.  7-1  are  shown  a 
sine  cycle  [ sin(2rrt/T)],  a cycle  of  the  same  energy  with  second 
harmonic  distortion  (.9  sin(2nt/T)  + .436  sin(4rrt/T))  , and  a 
cycle  of  the  same  energy  with  third  harmonic  distortion 
(.9  sin(2nt/T  + .436  sin(6nt/T)).  Consider  a sonar  pulse  made 
up  of  a set  of  sine  cycles  of  form  sin(2rrt/Tn)  . For  example,  for 
a Doppler  invariant  FM  slide,  Tfl  - TQ  + n • AT,  where  TQ  and  AT 
are  constants.  If  the  returning  signal  is  not  exactly  the  same 
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as  the  signal  sent  out,  a greater  processor  gain  will  be  obtained 

using  a replica  made  up  of  distorted  cycles  of  form  [a,  sin^^t/T  ) 

2 2 2 1 n 

+ a2  sin(4nt/Tn)  + a^  sin(6TTt/Tn)]  , ai  + a2  + a3  = 1,  than  for  a 

replica  of  the  same  as  the  outgoing  pulse.  The  requirement 
2 2 2 

al  + a2  + a2  * 1 will  insure  that,  for  a noise  input,  the  rms 
level  of  the  correlator  output  will  not  change  as  the  replica  is 
distorted.  Let  c^(t)  be  the  result  of  correlating  the  return 
with  the  outgoing  pulse,  c2(t)  be  the  result  of  correlating  the 
return  with  a pulse  made  up  of  second  harmonic  cycles  (sin  4nt/Tn) , 
and  c^t)  the  result  of  correlating  with  3rd  harmonic  cycles.  The 
correlation  of  the  return  with  the  distorted  replica  is  then 

3 3 

CD(t)  " Z ai  Ci(T)  ; X ai2  = 1 * (7.1) 

i=l  i^l 

The  problem  is  to  maximize  Cd(t)  by  proper  choice  of 

a^,  a2  and  a^.  A maximum  value  of  Cjj(t)  corresponds  to  a maximum 

value  of  Cd2(t) . Further,  it  is  noted  that  since  Cd(t)  and  Cd2(t) 

have  a rapid  fluctuation  as  r changes  it  would  be  of  more  signifi- 

cance  to  choose  the  a's  to  maximize  Cn  (t)  where  the  bar  denotes 

——k  u 

a time  average.  In  this  case  CD  (t)  corresponds  to  the  output  of 
a processor  that  consists  of  a linear  distorted  replica  correlator 
followed  by  a square  law  detector  and  a perfect  averager.  Let 
be  the  time  lag  that  the  processor  has  maximum  output  when  cor- 
relating the  output  with  the  input,  i.e.. 


time 


max . 


(7.2) 


2 

The  problem  is  then  to  choose  a^,  a2  and  a^  such  that  (T0)  = 

max.  Using  the  method  of  Lagrangian  multipliers  results  in 
three  equations. 
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al^Cl  " X)  + a2  ^1^2  + a3  = 0 

n 

a2  ^1^2  + a2^C2  ~ O + a3  ^2^3  = 0 

7 

a3  ^1^3  + a2  ^2^3  + a3^C3  " x)  = 0 • 


(7.3) 


Here  X is  the  Lagrangian  multiplier,  which  can  be  shown  to  equal 

CD  »(x  - Cp  ).  The  above  Eq.  (7.3)  has  the  trivial  solution, 

a^  - a£  = a^  * 0,  unless  the  determinant  of  the  coefficients  of 

the  a's  is  identically  zero.  Solution  of  this  determinant  leads 

to  three  values  of  X,  the  largest  of  which  corresponds  to  the 

maximum  correlator  output.  Substituting  this  value  of  X back 

into  any  two  of  the  equations  in  (7.3)  results  in  values  of  the 

2 2 2 

ratios  an<*  a3^ al*  Now  re<lu*-ring  that  a*  + + a^  = 1 

completes  the  solution 

The  utility  of  this  type  of  processing  depends  upon  the 
amount  of  improvement  and  whether  or  not  the  coefficients  a2  and 
a ^ measured  in  a large  number  of  cases  have  mean  values  signifi- 
cantly different  from  zero.  For  example,  if  for  1000  experimental 
pulses  maximum  correlation  was  obtained  for  a set  of  a2's  with 
mean  0.01  and  standard  deviation  0.1,  it  would  not  be  worthwhile 
to  correct  for  second  order  distortion  since  in  almost  half  of 
the  cases  setting  a2  * 0.01  would  reduce  the  processor  output. 
However,  if  the  mean  was  0.1  and  the  standard  deviation  was 
0.01,  then  the  use  of  this  correction  would  depend,  among  other 
factors,  upon  whether  or  not  the  processing  gain  was  significant. 

Distorted  replica  correlation  was  tried  on  7 samples 
of  the  5 msec  pulse  described  in  Section  5.  The  mean  value  of 
a2  was  calculated  to  be  ,004  with  a standard  deviation  equal  to 
.0003.  The  mean  value  of  a^  was  -.006  with  a standard  deviation 
equal  to  .001.  The  mean  processing  gain  for  the  "best"  replica 
in  each  case  was  only  .003  dB.  Because  these  calculations  are 
expensive,  in  terms  of  computer  time,  and  because  of  the  negligible 
gain  in  processor  output,  no  further  cases  were  run. 
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Although  the  result  of  this  effort  was  negative  the 
itsult  has  significance.  Modifying  the  shape  of  the  replica  in 
a significant  manner  had  essentially  no  effect  on  the  processor 
output.  Thus,  at  least  for  the  pulse  forms  and  ocean  surfaces 
tested  here,  there  was  no  significant  processing  gain  obtainable 
by  matching  the  return  to  a replica  with  second  and  third  order 
harmonic  distortion. 


34 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  78/21 


8.  IRREGULAR  SIGNAL  FORMS 

In  preceding  sections  sonar  pulses  have  been  formed 
by  stringing  together  a set  of  sine  cycles.  This  is  an  excellent 
representation  of  CW  or  FM  pulses  but  fails  to  represent  the 
pseudo- random  noise  pulse  or  an  explosive  pulse.  In  these  cases 
a generalization  of  the  preceding  theory  can  be  used. 

Consider  the  set  f^,  f2»  f^,  ...  obtained  by  sampling 
the  signal  f at  times  t^,  t2»  t^,  ...  where  t = nAt.  If  At  is 
chosen  such  that  the  sampling  rate  is  twice  the  highest  signifi- 
cant frequency  component  in  f then,  by  Shannon's  sampling  theorem, 
f is  completely  described  by  f^,  f2»  f ^ , ...  . The  next  step 
wi 11  be  to  fill  in  between  the  f^  with  a well  behaved  function 
that  is  amenable  to  integration  by  use  of  Kirchoff's  integral. 
Consider  the  points  f^,  f2  and  f^  shown  as  b's  in  Fig.  8-1.  A 
smooth  curve  of  form  A^  sin(2rrt/T^  + cp^)  + is  fit  between 
these  points.  A second  curve  of  the  same  form,  A0  sin(2nt/T2 
-4-  cp  2)  + B2 , is  fit  between  f^,  f^,  and  f^  with  the  restriction 
that  it  match  the  slope  of  curve  at  f^.  This  requires  solution 
of  a set  of  four  transcendental  equations  for  the  four  constants 
a2»  T2,  cp 2 and  b2*  It  may  be  that  the  solution  of  the  transcen- 
dental equations  can  only  be  made  in  the  least  mean  square  error 
sense.  If  the  least  error  is  not  acceptable  the  requirement  that 
the  curve  pass  through  f^  can  be  dropped  with  the  next  new  curve 
beginning  at  f^  rather  than  f^. 

The  form  of  the  above  curves  is  well  suited  for  use 
with  material  developed  in  earlier  sections  since  for  each  piece 
of  the  curve  [ f]  and  [ f']  have  simple  forms.  Thus  the  evaluation 
of  Eq.  (4.6)  is  not  much  more  complicated  than  it  is  for  f(t) 
made  up  of  a set  of  complete  sine  cycles. 

Other  mathematical  forms  can  be  used  to  represent  f(t) , 
for  example,  straight  line  segments.  Combinations  of  forms  are 
possible  such  as  an  explosive  pulse  represented  by  a linear 


FIG.  8-  1 - IRREGULAR  SIGNAL  FORM 
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pressure  rise  with  an  exponential  function  for  the  fall  off. 

Two  restrictions  must  be  kept  in  mind,  however.  First,  f(t) 
must  be  continuous.  Otherwise,  there  will  be  a delta  function 
contribution  from  the  [ f']  term  in  Eq.  (4.6).  Second,  when 
the  time  duration  of  a section  of  f(t)  is  decreased  there  may  be 
a corresponding  reduction  in  the  maximum  size  of  surface  elements, 
since  any  element  must  be  smaller  than  the  region  covered  by  any 
segment  of  f(t)  reflecting  from  the  surface. 
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9.  CONCLUSIONS 

This  report  covers  the  first  phase  of  a research  project 
to  study  the  correlation  loss  of  a sonar  pulse  when  it  reflects 
from  the  time-varying  ocean  surface.  A method  for  calculating  the 
reflected  pulse  has  been  developed,  based  on  Kirchoff's  integral, 
and  a computer  program  has  been  written  to  do  the  computations. 

Test  runs  on  a short  pulse  (5  msec)  show  an  average  correlation 
loss  of  1.6  dB  for  a pulse  reflecting  from  the  sea  driven  by  a 
10  knot  wind,  and  3.9  dB  for  a pulse  reflected  from  a 20  knot  sea. 
Almost  all  of  the  correlation  loss  is  due  to  an  energy  loss  in 
the  reflected  pulse. 

The  next  phase  of  this  work  will  concern  the  correlation 
loss  of  longer  pulses.  Other  interesting  and  useful  problems 
include  the  study  of  reflection  of  wide  beams,  the  correlation  loss 
of  pseudo-random  noise  signals,  the  addition  of  capillary  waves 
and  white  caps  to  the  sea  surface,  non-specular  scattering,  and 
so  on.  All  of  these  problems  can  be  studied  with  little  or  no 
modification  of  the  basic  theory  developed  here. 

Some  preliminary  calculations  for  100  msec  FM  slides  at 
1 kHz  and  3.5  kHz  have  been  run.  These  tend  to  show  slightly 
smaller  energy  and  correlation  losses  than  for  the  short  test 
pulses . 
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